Effect of a sub and suprathreshold periodic forcing an excitable glow discharge plasma 

near its bifurcation point. 
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In this paper nonlinear dynamics of a periodically forced excitable glow discharge plasma has been 
studied. The experiments were performed in glow discharge plasma where excitability was achieved 
for suitable discharge voltage and gas pressure. The plasma was first perturbed by a subthresh- 
old sawtooth periodic signal, and we obtained small subthreshold oscillations. These oscillations 
showed resonance when the frequency of the perturbation was around the characteristic frequency 
the plasma, and hence may be useful to estimate characteristic of an excitable system. On the 
other hand, when the perturbation was suprathreshold, system showed frequency entrainments. We 
obtained harmonic frequency entrainments for perturbation frequency greater than the character- 
istic frequency of the system and for lesser than the characteristic frequency, system showed only 
excitable behavior. 



I. INTRODUCTION 

Dynamics of a nonlinear plasmas under an external pe- 
riodic perturbations have been studied extensively in sev- 
eral plasma systems, where different nonlinear phenom- 
ena, like periodic pulling, frequency entrainment, etc., 
were observed [1-9]. For suitable plasma parameters 
these systems may also show excitable dynamics [9-11]. 
Excitable plasmas show noise-induced resonances under 
both stochastic and periodic perturbations [9-11] that 
have also been observed in many other physical, chemi- 
cal, biological and electronics systems [12-27]. Though 
excitable systems show rich dynamical behavior under a 
pure periodic forcing, little study is carried out in this 
area [28, 29]. In this paper, our aim to explore the dy- 
namics of glow discharge plasma near excitable region 
due to periodic perturbations. 

The basic characteristic of an excitable system is that 
it shows a fixed point or coherent limit cycle oscillations 
depending on the value of the control parameter (CP) of 
the dynamics. The point where the system changes from 
oscillatory to fixed point behavior is called threshold or 
bifurcation point. Now if the system shows fixed point 
behavior for the value of CP below the threshold, then 
the dynamics will be limit cycle oscillations on the other 
side of the threshold or vice versa, depending upon the 
system properties. Response of an excitable system at 
fixed point state to an external perturbation applied on 
the CP depends on the perturbation amplitude. When 
the perturbation is subthreshold, i.e., the amplitude is 



* jaman_nonlinear@yahoo.co.in 
t ansekar.iyengar@saha.ac. in 



small such that the threshold is not crossed, the sys- 
tem shows small oscillations around the fixed point [23]. 
These subthreshold or small oscillations may bear the 
characteristic of system dynamics [30-32]. On the other 
hand, when the perturbation is suprathreshold, i.e., it's 
amplitude is large enough to cross the threshold, the sys- 
tem returns to its fixed point deterministically, i.e., once 
the threshold is crossed, the system becomes almost inde- 
pendent of the perturbation and comes to its fixed point 
state traversing one limit cycle [13-16, 33]. In plasma, 
excitability may be obtained through Hopf [9, 11] or ho- 
moclinic bifurcation [10] depending on the plasma pa- 
rameters. 



Usually, dynamics of a periodically perturbed plasma 
depends on both the amplitude and frequency of the ex- 
ternal perturbation that may be represented by "Arnold 
tongue" diagram. Though the dynamics of an excitable 
system depends on the frequency of the perturbation, it 
does not depend much on the perturbation amplitudes 
as long as the perturbation lifts the system at its excited 
state. This is because, the variation in perturbation am- 
plitude does not affect much the system dynamics at ex- 
cited state [23]. Hence usual "Arnold tongue" diagram 
may not be obtained. Here, we have presented the dy- 
namics of an excitable glow discharge plasma under sub- 
threshold and suprathreshold periodic sawtooth forcing. 
The paper has been organized as follows: we present a 
brief description of the experimental setup in Sec. II and 
autonomous system dynamics through which excitabil- 
ity was obtained in the glow discharge plasma has been 
discussed in Sec. III. Experimental results has been pre- 
sented and discussed in Sec. IV. Finally, we conclude the 
results in Sec. V. 
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II. EXPERIMENTAL SETUP 
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FIG. 1. Schematic diagram of the cylindrical electrode sys- 
tem of the glow discharge plasma. The probe is placed at a 
distance / ~ 12.5 mm from the anode. 

The experiments were carried out in a glow discharge 
argon plasma produced by a dc discharge in a cylindrical 
hollow cathode electrode system with a typical density 
and temperature 10 7 cto -3 and 2 — 6 eV respectively 
The schematic diagram of the electrode system is shown 
in Fig. 1. The electrode assembly was housed inside 
a vacuum chamber and the neutral pressure inside the 
chamber was controlled by a needle valve and the range of 
the gas pressure in these experiments was between 0.001 
to a few mbar. The discharge voltage (0 — 999 V) was 
applied between the cathode and the anode keeping the 
anode grounded. A signal generator (Fluke PM5138A) 
was coupled with the power supply for the experiments of 
periodically forced plasma. The main diagnostics was the 
Langmuir probe which was used to obtain the floating po- 
tential fluctuations. In these experiments the discharge 
voltage (DV) and pressure were the control parameters. 
Detail of the experiments can be found in Ref. [10]. 



the glow unstable [34] , that was probably responsible for 
the appearance of these oscillations [7, 35] . 

Increasing the DVs, lead to an increase both in the 
amplitude and the time period of the oscillations and 
they become regular (period 1) at certain point of DV. 
Further increase in the DV modifies the oscillation profile 
and results in typical relaxation oscillations which change 
to excitable fixed point behavior at certain point. The 
DV at which these oscillations cease may be termed the 
bifurcation point (Vh) which depends also on the neutral 
pressure. The higher the neutral pressure, the higher the 
Vh- As the frequencies of the instabilities were within 
the ion plasma frequency (order of a few kHz) and the 
presence of anode glow and relaxation oscillations have 
been attributed to the presence of double layers, these 
instabilities could be related to ion acoustic instabilities 
and also anode glow related double layer [36, 37]. 

The time period (T) of these relaxation oscillations in- 
creases with increasing DV near the bifurcation point, 
and becomes infinite beyond the Vh, results in vanish- 
ing of the limit cycles. For larger values of DV, the au- 
tonomous dynamics exhibits a steady state fixed point 
behavior. Time traces from top to bottom in the left 
panel of Fig 2 depict this period lengthening of the oscil- 
latory behavior. A systematic analysis of the increment 
in the period (T), presented in Fig. 2(a) [right panel], 
indicates that the autonomous dynamics undergoes a ex- 
ponential slowing down. Consequently, the ln|V — Vh\ 
vs T curve can be fitted by a straight line. The results 
of Fig 2 indicate that the system dynamics undergo a 
homoclinic bifurcation at Vh [10]. The point Vh divides 
the plasma into two distinct regions. The floating poten- 
tial fluctuation exhibits relaxation oscillations on the one 
side of the Vh and stable fixed point on the other side, 
which is a typical behavior of an excitable system. In the 
experiments, DV was so chosen that the system exhibited 
fixed point behavior and were perturbed by subthreshold 
and suprathrcshold external sawtooth periodic signals. 

IV. RESULTS 

Dynamical behavior of the system under such pertur- 
bation has been presented in this section. 



III. AUTONOMOUS DYNAMICS 

By varying the neutral pressure, the breakdown of the 
gas occurs at different DV. At higher pressure (i.e., or- 
der of 1 mbar), with the initiation of discharge, a bright 
glow or spot was observed to form on the anode that was 
unstable initially, and was the source of the fluctuations. 
An interesting feature associated with the anode glow 
was the different types of oscillations in the floating po- 
tential at different pressures. At the initial stage, system 
showed irregular and relaxation type oscillations. Ions 
produced inside the anode glow due to collisions of the 
accelerated electrons across it with the neutrals, made 



A. Effect of subthreshold signal 

For the study of the effect of subthreshold signal, the 
reference voltage Vq was chosen such that Vq > Vh and 
therefore the autonomous dynamics, by virtue of an un- 
derlying homoclinic bifurcation, exhibits steady state be- 
havior. The discharge voltage V was thereafter per- 
turbed, V = Vq + S(t), where, the subthreshold saw- 
tooth periodic signal S(t) was so chosen that the volt- 
age V = Vq + S(t) < Vh, i.e., the applied perturbation 
does not cause the system to cross over to the oscillatory 
regime. The subthreshold periodic sawtooth signal of 



3 



0) 

3 
o 

> 



10 



-10 



\\\\\\\\\\\n 



; k iv i,ii> | 

u Ika | liwkniim ^WtalL 

npp '^W ™PI 



0.20 


(a) 


I 


_ °- 15 

£ o-io 




/ 
p 


H- 1 0.05 




> 


0.00 




6 


300 350 


400 450 500 550 6 
DV(V) 


5 
4 


V 

■ 


1 i | i | 

(b) " 








£ 2 

t 







0.05 0.10 0.15 0.20 0.25 

time (ms) 



0.5 1 

Time (ms) 

FIG. 2. Left panel: Time series showing that the period of relaxation oscillations increases with augmenting DV. Right panel 
:(a) Exponential increment of the time period (T) with DV and (b) In | V — Vh\ vs T curve can be fitted by a straight line 
indicating an underlying homoclinic bifurcation. 



frequency, and showed maximum value at frequency 3.2 
kHz. Power again decreased with increase in perturba- 
tion frequency, hence system showed resonance at 3.2 
kHz. As the frequency of the internal plasma oscilla- 
tions which were basically the acoustic instabilities, was 
of the order of a few kHz in these parametric regions of 
plasma [10, 36, 37], so system shows resonance at inter- 
nal frequency. Therefore, these subthreshold oscillations 
may be useful to determine the system frequency. One 
important observation was that these oscillations were 
not observed far away from the bifurcation point. This 
indicates that the oscillations are the feature of an ex- 
citable dynamics [23]. This approach may be useful to 
determine the frequency of internal oscillations of other 
excitable systems. 
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FIG. 3. Resonance curve of the subthreshold oscillations with 
frequency of the applied sawtooth signals. It peaks at system 
internal frequency around 3.2 kHz. 



variable frequency was generated using a Fluke PM5138A 
function generator. In this experiment, the frequency of 
the applied signal was varied from a few Hz to 20 kHz. 
When the frequency of the applied signal was around a 
few kHz, system showed subthreshold oscillations. Am- 
plitude of these oscillations was increased with increase 
in the perturbation frequency that was recorded using an 
oscilloscope. 

Figure 3 shows the power of the subthreshold oscilla- 
tions as a function of the frequency of the applied signals. 
It shows that initially the power of the subthreshold os- 
cillations was minimum and increased with increase in 



B. Effect of suprathreshold signal 

For the experiments on the effect of a suprathreshold 
forcing, the DV (Vq) was set such that the floating po- 
tential fluctuations exhibit fixed point behavior. In order 
to minimize the effect of parameter drift, a set point (Vq) 
quite far from the homoclinic bifurcation (Vh) was cho- 
sen. Subsequently, a suprathreshold sawtooth signal of 
fixed amplitude but variable frequency was superimposed 
on the DV. So in this case perturbation always crosses 
the threshold, i.e., V = Vo+S(t) > Vh- The frequency of 
the applied signals was varied from a few Hz to 20 kHz. 

At low frequency, i.e., less than 1 kHz, though the 
amplitude of the applied signal is sufficient to cross the 
threshold, the system did not show any excitable dynam- 
ics. In these frequency, we obtained only the plasma 
distorted sawtooth signals as shown in Fig. 4(a) for per- 
turbation frequency 0.2 kHz. We obtained perturbation 
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FIG. 4. Left panel: time series of plasma floating potential for periodic sawtooth perturbation: (a) at frequency 0.2 kHz; (b) 
at 1.072; and c) at 3.56 kHz, respectively, and their power spectra at (a'), (b'), and (c') respectively. Figure(a) shows that 
system does not responds at frequency 0.2 kHz; figure(b) shows that the system its excitable behavior at 1.072 kHz and Fig.(c) 
shows 1:1 entrained state at 3.56 kHz. 1:1 entrained state is also clear from the power spectrum [Fig.(c')] also where internal 
frequency of the system coincides with external perturbation. 



invoked limit cycle oscillations at frequency 1.072 kHz 
as shown in Fig. 4(b). In this figure, the oscillations 
with larger and smaller time periods are the applied saw- 
tooth signal and the limit cycle oscillations of the plasma 
system at excited state respectively. Oscillations with 
smaller time period were of the order of system's charac- 
teristic time scale (inverse of internal frequency). From 
Figs 4(a) and (b), it is also clear that the excitation mech- 
anism of an excitable system depends not only on the 
perturbation amplitudes but also on the frequency. 

Fig. 4(c) and (c') show the system output and cor- 
responding power spectrum when the perturbation fre- 
quency was 3.56 kHz. They clearly show that the system 
response is phased locked with the external periodic sig- 
nal with rotation number 1:1. In case of phase locked 
states, rotation number (w) is defined as W = f/fo, 
where / and fo are the frequencies of external pertur- 
bation and system frequency, respectively [8]. So when 
the oscillations frequency of the system was close to per- 
turbation frequency (/ ~ /o) that is also clear from the 
frequency measurements presented in Subsec. IV A, 1:1 
phase locking or entrainment phenomenon was obtained. 
Figures 4(b) and (c) show that the periodic perturbations 
invoked limit cycle oscillations at frequencies 1.072 and 
3.56 kHz respectively. The results shown Figures 4(b) 
cannot be taken as frequency entrainment phenomenon, 
as the phase is not locked with the external perturba- 



tion. In this case, external perturbation lifts the system 
only at excited state and hence shows limit cycle oscilla- 
tions. Whereas, at 3.56 kHz [Fig. 4(c)], external pertur- 
bations lifts the system at excited state as well as locks 
the phase of internal oscillations. From these observa- 
tions it is clear that behavior of excitable plasma under 
external perturbation depends both on its excitable and 
internal dynamics. 

When the perturbation frequency was increased to 8.93 
kHz, system showed 1:2 harmonic frequency entrainment 
as shown in Fig. 5(a) and corresponding power spectrum 
shown in Fig. 5(a'). In Fig. 5(a'), fo represents the exter- 
nal perturbation and fo/2 is the system oscillation at ex- 
cited state. At perturbation frequencies 9.48, 14.38, and 
20.0 kHz, we obtained 1:3, 1:4, and 1:5 entrainments or 
phase locked states, respectively. The entrainment phe- 
nomena is also clear from the times series shown Figs. 5 
(b), (c) and (d). These facts are also clear from their 
power spectra shown in Fig 5 (b'), (c') and (d') respec- 
tively. 

Between any two phase locked states, system re- 
sponded such that it showed intermittent locked states. 
Figures 6(a) and (b) show the typical intermittent pat- 
terns between 1:2 and 1:3; and 1:4 and 1:5 respectively. 

The above results show that the low frequency per- 
turbations (less than its intrinsic frequency) only lift the 
system at excited state only, and hence system shows 
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FIG. 5. Time series of plasma floating potential for periodic sawtooth perturbation: (a) at frequency 8.93 kHz; (b) at 9.48; c) 
at 14.38 kHz; and d) at 20.0 kHz respectively. Figure(a) shows 1:2 entrained; figure(b) shows 1:3 entrained state; Fig.(c) shows 
1:4 entrained state; and Fig.(d) shows 1:5 entrained state. 



V. CONCLUSION 
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In this paper we have studied the dynamical behav- 
ior of excitable glow discharge plasma under the sub- 
threshold and suprathrcshold periodic forcing. System 
shows subthreshold oscillations under subthreshold forc- 
ing. These oscillations are the characteristics of the sys- 
tem dynamics and show resonance with external pertur- 
bations at the internal frequency of the system. Hence 
such resonance curve may be useful to characterize ex- 
citable plasma and this idea may also be applied in other 
excitable systems. 



limit cycle oscillations which are not phase locked with 
perturbations. But for plasmas which exhibit van der Pol 
type behavior, show 2:1, 3:1, etc., subharmonic entrain- 
ments under low frequency [2, 8]. In excitable plasma 
driving amplitude does not play not much role in the 
dynamics as long as the amplitude is sufficient to excite 
the system and hence Arnold tongue diagram may not 
be obtained [2, 8]. 



For the suprathrcshold perturbation, system showed 
excitable or mode locking behavior depending on the 
perturbation frequency. Most simple frequency entrain- 
ment (1:1) phenomenon was obtained when the driver 
frequency was chosen close to system frequency. When 
the driver frequency was less that the system frequency, 
system showed natural excitable dynamics, whereas, for 
higher frequency, harmonic entrainments were observed. 
As the dynamics of an excitable system at excited state 
is almost independent of the amplitude of external per- 
turbation, we did not get Arnold tongue behavior in the 
excitable region of glow discharge plasma. 
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